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Abstract —  A  radar  pulse  that  impinges  upon  a  radar  resolution 
area  of  the  sea  snrface  produces  backscattered  returns  which  are 
called  sea  clutter.  For  low  grazing  angles  and  very  fine  radar 
cell  resolution  areas,  the  clutter  intensity  distribution  departs 
significantly  from  the  exponential  distribution.  In  this  case,  the 
clutter  is  said  to  display  spiky  behavior  and  the  distribution  of  the 
intensity  develops  a  much  longer  tail  relative  to  the  exponential 
distribution.  Statistical  analysis  of  collected  data  near  a  grazing 
angle  of  0.2  °  at  X-band  from  the  sea  off  the  coast  of  Kauai, 
Hawaii  are  examined  relative  to  the  log-normal,  Weibull  and 
K  distributions.  Based  on  an  analogy  of  sea  clutter  and  other 
disciplines  including  computer  networks  and  finance,  we  also 
apply  the  Pareto  distribution  to  the  collected  data.  We  also 
compare  the  data  to  the  WW  and  KK  two-population  mixture 
distributions.  Maximum  likelihood  estimation  of  the  parameters 
of  the  dlstrlbntions  are  obtained  from  the  measured  data.  In 
all  cases,  the  two  population  mixture  distributions  and  the 
Pareto  distribution  are  more  accnrate  than  the  three  classical 
distributions.  However  the  Pareto  distribution  has  the  advantage 
of  being  an  analytically  tractable  two  parameter  distribution 
while  having  similar  accuracy  to  the  five  parameter  WW  and 
KK  at  critical  valnes. 

1.  Introduction 

The  performance  of  high  resolution  radars  at  low  grazing 
angles  has  been  a  topic  of  research  for  the  last  several  decades. 
Extensive  surveys  of  many  of  the  experimental  and  theoretical 
developments  are  given  in[l][2].  Despite  the  extensive  number 
of  studies,  the  current  phenomenological  understanding  of  sea 
clutter  modeling  is  incomplete.  Partially  due  to  this  incomplete 
knowledge,  a  variety  of  statistical  models  have  been  proposed 
for  the  distribution  of  the  sea  clutter.  An  accurate  sea  clutter 
model  allows  for  the  optimization  of  the  radar  detector  chain 
including  its  GEAR  processor.  This  in  turn  allows  for  enhanced 
target  detection  and  better  control  of  the  false  alarm  rate. 
In  this  paper  we  assess  the  accuracy  of  some  new  and 
unexamined  mixture  distributions  for  the  low  grazing  angle  sea 
clutter  data  of  [4].  Maximum  likelihood  estimation  (MLE)  of 
the  parameters  of  the  distributions  are  obtained  and  the  results 
are  compared  to  the  actual  data’s  distribution.  In  Section  II 
of  this  paper  we  briefly  review  the  measurement  parameters, 
the  collected  data  and  sea  conditions  of  the  experiment.  In 
Section  III,  we  overview  the  traditional  models  as  well  as  the 
mixture  models[5]  and  Pareto  distribution [7].  In  Section  IV 
we  estimate  the  parameters  of  the  distributions  using  the  MLE 
and  compare  their  respective  CDEs.  Conclusion  of  this  paper 
is  given  in  Section  V. 
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11.  Experimental  Measurements 

An  itemized  description  of  the  measurement  radar  appears 
in  Table  1.  The  Data  was  collected  in  Kauai,  Hawaii  between 
August  and  December  of  1994  by  the  Naval  Research  Labo¬ 
ratory  (NRL),  in  conjunction  with  Naval  Air  Warfare  Center, 
China  Lake  (NAWC/CL),  the  Johns  Hopkins  University  Ap¬ 
plied  Physics  Laboratory  (JHU/APL),  and  Texas  Instruments 
(TI),  since  then  acquired  by  Raytheon.  There  are  several 
descriptors  to  note  about  the  data.  The  grazing  angle,  which 
is  a  function  of  the  radar  height  and  radar  geometry  is  around 
0.22  degrees.  The  data  was  collected  for  both  VV  and  HH. 
However  the  Radar  was  not  dual-polarized.  Consequently,  the 
HH  and  VV  data  collections  were  measured  separately.  The 
Radar  Height  was  23  meters  and  the  radar  range  was  5.74 
km  for  the  vertical/vertical  (VV)  measurements  and  6.11  km 
for  the  horizontal/horizontal  (HH)  measurements.  The  second 
major  descriptor  of  the  measured  data  is  the  transmit  geometry. 
Both  upwind  (UP)  and  crosswind  (CR)  measurements  were 
obtained.  However,  in  this  paper,  we  have  only  concentrated 
on  the  upwind  (UP)  measurements  as  they  exhibit  more  of 
the  ’’spiky”  characteristics  which  makes  sea  clutter  difficult 
to  model  at  low  grazing  angles.  During  the  course  of  the 
clutter  measurements,  the  wind  speed  measured  was  about  9 
m/s,  while  the  largest  significant  wave  height  (average  height 
of  the  1/3  highest  waves)  measured  was  around  3  m.  These 
conditions  indicate  a  hydrographic  sea  state  of  4.  In  Eigure 
1,  we  have  presented  the  calibrated  RCS  in  decibels  above 
a  square  meter  (dBsm)  as  a  function  of  time  in  seconds 
along  the  horizontal  axis  and  of  range,  in  meters,  along  the 
vertical  axis.  These  are  also  known  as  range-time  intensity 
(RTI)  plots.  In  Eigure  2,  we  have  plotted  the  variation  of 
the  HH  and  VV  RCS  for  a  specific  range  cell  over  a  5 
seconds.  As  noted  in  [4],  the  HH  clutter  is  sharper,  spikier 
and  more  intermittent.  In  Eigure  3,  the  temporal  variation  of 
the  spectral  frequency  content  is  displayed  for  the  range  cell 
of  Eigure  2  in  the  following  manner.  The  individual  spectra  are 
generated  by  diving  the  entire  time  interval  into  consecutive, 
non-overlapping  windows,  and  then  calculating  the  power 
spectral  density  (PSD)  within  each  Hanning-weighted  window. 
The  set  of  power  spectra  is  then  normalized  by  dividing  each 
value  in  every  PSD  by  the  largest  value  that  occurred  in  entire 
set;  the  new  maximum  is  1.  Einally,  the  entire  sequence  of 
normalized  power  spectral  spanning  the  entire  time  interval  is 
displayed  in  a  3-D  plot. 
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vertically  polarized,  upwind  sea  clutter 


Fig.  1.  RTI  plots  of  radar  cross  section  in  dB  above  a  square  meter  for  VV 
and  HH 
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Fig.  2.  Radar  cross  section  in  square  meters  for  a  specific  range  of  the  HH 
and  VV  data 
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vertically  polarized,  upwind  sea  clutter:  example  one.  observed  with  50  millisecond  window 
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Fig.  3.  Normalized  power  spectra  for  HH  and  VV  for  a  specific  range  cell 


TABLE  I 

Radar  Descriptors 


Operation  Frequency 

X-band  (9.5-10.0  Ghz) 

Average  Power 

500  W 

Pulse  Length 

2.5  ns  compressed 

Pulse  Repitition  Frequnecy  (PRF) 

2000  Hz 

Data  Collection  Extent 

156  m  range  swath 

Data  Collection  Mode 

Fixed  antenna  pointing 

Azimuth  Beam  width 

2.4° 

Elevation  Beam  width 

4.0° 

Polarization 

VV  and  HH 

III.  Distribution  models 

A  variety  of  distribution  models  have  been  proposed  to 
model  distribution  of  the  sea  clutter  amplitude  or  intensity 
returns.  Among  these,  the  K[l],  Log-Normal[3]  and  Weibull[2] 
distributions  have  been  cited  frequently  and  partial  success 
with  these  models  have  been  reported.  From  these,  the  K 
distribution  which  is  a  compound-Gaussian  model)  13]  [14] 
has  been  given  a  phenomenological  interpretation.  According 
to  this  interpretation,  the  K  distribution  envelope  is  a  com¬ 
pound  distribution  consisting  of  a  locally  Rayleigh  distribution 
speckle  whose  mean  is  modulated  by  a  gamma  distribution. 
The  speckle  part  is  assumed  to  be  due  to  the  scattering 
primarily  by  the  capillary  waves  and  ripples  while  the  under¬ 
lying  mean  (also  called  the  texture)  of  the  patch  is  assumed 
to  be  from  the  ocean  gravity  waves.  Instead  of  the  gamma 
distribution,  other  distributions  for  the  texture  such  as  log¬ 
normal,  generalized  Gaussian  and  inverse  gamma  have  been 
proposed  [13]  [14].  In  their  work  [13],  it  was  noted  that  for  a 
range  resolution  of  3  meters,  the  compound-Gaussian  and  the 
K  distribution  did  not  provide  a  good  fit.  The  K  distribution 
has  also  been  extended  to  take  into  account  different  types 
of  ’’spiky”  phenomenon)  1]  such  as  discrete  spikes  (bursts 
and  whitecaps).  For  example,  the  compound  K-A  distribution 
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consists  of  the  Poisson  distribution,  gamma  distribution  and 
the  Rayleigh  distributions[l].  Taking  a  different  approach, 
Dong  has  recently  [5]  proposed  the  KK  and  Weibull-Weibull 
(WW)  distributions  for  high  resolution  and  high  grazing  angle 
sea  clutter.  His  results  indicates  that  the  KK,  WW  and  KA 
perform  similarly  in  such  a  scenario.  However  to  estimate 
the  parameters  of  the  KK  and  WW  is  simpler  than  the 
K-A  distribution.  These  mixtures  distribution  can  be  given 
the  interpretation  of  a  two-population  model.  That  is  with  a 
probability  of  p  the  clutter  belongs  to  a  certain  distribution 
with  a  CDF  Fi{z)  and  with  a  probability  of  1  —  p,  the  clutter 
belongs  to  another  distribution  with  a  CDF  F2(z).  Here,  we 
briefly  overview  the  distributions  used  in  this  paper. 

The  Probability  Density  Function  (PDF)  of  the  log-normal 
distribution  is  given  by[3] 


fUz) 


1 


- =  exp- 

zcrv27r 


{Injz)  -  pf 
2a2 


(1) 


where  z  is  lognormal  if  and  only  if  ln{z)  is  a  Gaussian  with  a 
mean  p,  and  variance  cr.  The  Cumulative  Distribution  Function 
(CDF)  FL{a)  is  of  the  form: 


F,(a)  =  m'ZFIzJL) 


(2) 


where  is  the  standard  normal  function. 

The  PDF  of  the  K  distribution  is  given  by[l]: 


fK{z)  = 


07Z±i 
26  2  2:  2 


7(i 


(3) 


where  z  denotes  the  intensity  of  the  sea  clutter,  v  is  the  shape 
parameter,  b  is  the  scale  parameter.  The  CDF  Fk  (a)  is  of  the 
form: 


Fk(z)  =  1  -  ^  Ky{2Vbz) 

li’Z) 


(4) 


Details  about  the  K  distribution  can  be  found  in  [1]. 

The  Weibull  distribution  which  contains  the  Rayleigh  dis¬ 
tribution  as  a  special  case  is  another  popular  distribution  for 
sea  clutter[3].  The  PDF  of  the  Weibull  distribution  is  given 
by[l]: 

fw{z)  =  exp{-{—y)  (5) 

UJ  'LU  'UJ 

where  z  denotes  the  intensity  of  the  sea  clutter,  7  is  the  shape 
parameter,  vu  is  the  scale  parameter.  The  CDF  Fw{z)  of  the 
Weibull  is: 


Fw{z)  =  l-exp{-{—y)  (6) 

zu 

The  WW  distribution  has  the  PDF: 


The  KK  distribution  is  defined  similarly,  except  the  mixtures 
are  K  distributions  instead  of  Weibull  distributions. 

The  PDF  of  the  Generalized  Pareto  distribution  is  dehned 
as[7]: 

/gp(^)  =  ^(1  -  (9) 

where  k  is  the  shape  parameter  and  A  is  the  scale  parameter. 
The  value  for  z  is  0  <  z  <  00  for  k  <  0  and  0  <  z  <  I- 
for  A:  >  0.  When  k  =  0,  the  Generalized  Pareto  distribution 
takes  the  form  of  an  exponential  PDF  with  the  mean  A.  When 
A:  <  0,  the  generalized  Pareto  distribution  takes  the  form  of 
Pareto  distribution.  When  A:  =  1,  it  takes  the  form  of  a  uniform 
random  variable  on  [0,  A].  For  this  paper,  k  is  always  less  than 
zero  and  z  denotes  the  intensity.  The  Pareto  distribution  is  also 
a  continuous  mixture  of  exponential  distributions  with  gamma 
mixture  parameters  and  it  is  formulated  as: 


fp{z;k)=  fE{z;OfG{yk,X)d^  (10) 

Jo 

where  Je  is  an  exponential  PDF  with  a  mean  ^  and  fc  is 
a  gamma  pdf  with  parameters  k,  A.  From  a  derivation  point 
of  view,  the  difference  between  the  Pareto  distribution  and 
the  K  distribution  is  that  the  reciprocal  of  the  mean  of  the 
exponential  in  the  Pareto  distribution  is  conditioned  by  the 
gamma  distribution.  In  the  K  distribution,  it  is  the  mean  of 
the  Rayleigh  (for  envelope)  and  exponential  (for  intensity) 
that  is  conditioned  by  the  gamma  [1].  The  Pareto  distribution 
has  some  shared  properties  including  the  power-law  property 
of  the  Tsallis  distribution  which  was  proposed  in  [12]  for 
modeling  the  differentiated  data  of  sea  clutter. 

It  should  be  noted  that  the  Pareto  distribution  has  only  two 
parameters  compared  to  the  five  parameters  of  the  WW  and 
KK  distributions.  The  Pareto  distribution  has  been  used  for 
heavy  tailed  problems  in  different  fields  including  physics, 
economics  and  finance[ll].  In  computer  networks,  it  has 
been  shown  that  data  traffic  may  exhibit  bursts  of  high  data 
rates  even  when  the  average  data  rate  is  low[10].  The  Pareto 
distribution  is  more  accurate  than  the  traditional  Poisson  and 
Binomial  distributions  for  such  a  case  [10].  This  is  analogous 
to  sea  clutter  where  the  average  clutter  return  could  be  low 
but  bursts  of  spikes  are  exhibited.  An  interpretation  is  that 
the  Pareto  distribution  follows  the  power  law[ll],  where  the 
probability  of  measuring  a  value  varies  inversely  as  a  power 
(with  an  exponent  factor)  of  that  value. 

IV.  Analysis  of  the  Data 


fww{z)  =  pfwi{z)  +  (1  -  p)fw2iz)  (7) 

and  its  CDF: 

Fww{z)  =  pFwiiz)  +  {1  -  p)Fw2{z)  (8) 

where  the  subscripts  indicate  two  Weibull  distributions  that 
may  have  different  shape  and  scale  parameters. 


To  obtain  the  parameters  for  a  distribution  fz  from  the  data, 
we  have  used  maximum  likelihood  estimator  (MLE)[8].  In 
[14],  it  was  shown  that  for  realistic  sea  clutter,  the  MLE  has  a 
smaller  mean  square  error  (MSE)  than  the  method  of  moments 
(MoM)  and  method  of  fractional  moments  (MoEM)  estimators. 

Eor  the  case  of  independent  and  identically  distributed 
samples,  if  N  independent  samples  zi,  Z2...,  z„  are  drawn  from 
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a  distribution  with  parameters  6,  then  their  likelihood  function 
is  defined  as: 

N 

L{0)  =  Y\_  fz{Zij0ML)  (11) 

n—1 

The  maximum  likelihood  estimate  of  8  denoted  by  6m l  is 
given  by: 

^ML  =  argmaxL(0)  (12) 

0 

In  order  to  ensure  that  the  joint  PDF  can  be  factored  into 
independent  and  identically  distributed  PDFs,  the  data  was 
sub-sampled  uniformly  by  a  factor  much  greater  than  its 
significant  correlation  length.  However,  to  ensure  accuracy, 
the  number  of  data  points  used  for  the  estimation  was  near 
a  half  million  samples.  For  the  Log-Normal  distribution,  the 
MLE  has  a  closed  form  solution.  For  the  K,  WW  and  KK,  we 
have  used  the  constrained  optimization  method  noted  in  [6].  In 
all  cases,  the  Iog-likelihood[8]  which  gives  numerically  more 
stable  results  were  used.  For  the  Generalized  Pareto  method, 
the  MLE  may  be  found  by  a  constrained  optimization  problem 
using  the  Nelder-Mead  aIgorithm[9]. 

Before  calculating  the  MLE  of  the  parameters,  we  normal¬ 
ized  the  clutter  intensity  by  one.  This  does  not  change  the 
shape  parameter  of  the  estimators  but  was  mainly  done  so  that 
the  slope  of  the  HH  and  VV  CDFs  can  be  compared  on  the 
same  axis.  Eor  the  actual  RCS  values  of  the  data,  the  reader 
is  referred  to  [4].  Figures  4  and  5  shows  a  plot  of  the  CDFs 
with  the  MLE  parameters  for  the  HH  data.  Eigures  6  and  7 
shows  the  CDEs  with  the  MLE  for  the  VV  data. 

In  Tables  2  and  3  we  have  calculated  the  difference  in  dB 
between  the  true  and  estimated  intensity  for  some  intermediate 
and  critical  CDE  values.  Eor  the  case  of  two  parameter 
distributions,  the  Pareto  has  the  best  fit  for  the  HH  and  VV 
data.  It  also  has  the  smallest  maximum  deviation  among  the 
two  parameter  distributions.  Compared  to  the  WW  and  KK 
distributions,  at  a  CDE  value  of  0.99,  the  Pareto  is  about  3 
dB  worse  for  HH  and  1.5  dB  worse  for  VV.  However  at  CDF 
value  of  0.999,  it  performs  better  than  the  WW  and  has  similar 
performance  to  the  KK.  We  should  note  that  WW  and  KK 
are  five  parameters  distributions  whereas  the  Pareto  is  a  two 
parameter  distribution.  Consequently,  one  would  expect  that 
a  distribution  with  five  parameters  outperforms  a  distribution 
with  two  parameters.  However,  the  Pareto  distribution  outper¬ 
formed  all  the  two  parameter  distributions  and  its  performance 
is  virtually  the  same  as  the  WW  and  KK  at  the  critical  point 
of  0.999. 

V.  Conclusion 

Several  distributions  were  examined  for  high  resolution 
and  low  grazing  angle  sea  clutter  data.  MLE  estimation  was 
performed  to  obtain  the  parameters  of  all  the  distributions. 
Despite  the  fact  that  the  Pareto  distribution  is  a  two  parameter 
distribution,  it  outperformed  the  W,  K,  log-normal  and  WW 
distributions  in  the  tail  region  and  its  performance  was  very 
close  to  the  KK.  We  also  note  that  the  Pareto  distribution 


TABLE  II 

CDF  Estimation  errors  for  HH  in  dB 


CDF  value 

0.25 

0.5 

0.75 

0.90 

0.99 

0.999 

Weibull 

2.02 

1.08 

2.88 

2.23 

4.72 

8.72 

Lognormal 

0.20 

0.61 

1.17 

0.30 

5.63 

8.12 

K 

0.21 

1.54 

2.55 

1.5 

5.74 

9.79 

WW 

0.39 

0.30 

0.14 

0.21 

0.52 

3.31 

KK 

0.21 

0.21 

0.19 

0.28 

0.37 

2.81 

P 

0.65 

0.21 

1.00 

0.53 

3.39 

2.74 

TABLE  III 

CDF  Estimation  ERRORS  for  VV  in  dB 


CDF  value 

0.25 

0.5 

0.75 

0.90 

0.99 

0.999 

Weibull 

1.45 

0.71 

1.78 

1.38 

3.15 

7.57 

Lognormal 

0.17 

0.07 

0.54 

0.37 

2.57 

5.16 

K 

0.45 

0.77 

0.86 

1.4 

3.60 

7.83 

WW 

0.17 

0.23 

0.06 

0.27 

0.14 

3.10 

KK 

0.17 

0.12 

0.11 

0.20 

0.21 

2.50 

P 

0.59 

0.05 

0.52 

0.37 

1.73 

2.53 

CDF  for  Horizontal  Polarization 


Fig.  4.  CDF  plots  of  the  actual  data  and  estimated  Weibull,  K  and  lognonnal 
distributions  HH 


CDF  for  Horizontal  Polarization 


Fig.  5.  CDF  plots  of  the  actual  data  and  estimated  WW,  KK  and  pareto 
distributions  for  HH 
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CDF  for  Vertical  Polarization 


Fig.  6.  CDF  plots  of  the  actual  data  and  estimated  Weibull,  K  and  lognonnal 
distributions  VV 
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Fig.  7.  CDF  plots  of  the  actual  data  and  estimated  WW,  KK  and  pareto 
distributions  for  VV 


is  analytically  tractable  and  simpler  than  the  K  distribution 
as  well  as  the  compound-Gaussian  distributions  in  [13]  [14]. 
For  the  analysis  of  some  radar  systems,  this  could  possibly 
lead  to  closed  formed  solutions  further  down  in  the  detection 
chain.  The  Pareto  distribution  has  also  been  noted  in  computer 
networks  and  other  fields  for  a  phenomenon  that  is  analogous 
to  sea  clutter.  While  a  general  trend  has  been  pushing  sea 
clutter  into  multiple  parameter  distributions  (five  in  the  case 
of  WW,  KK  and  KA),  the  two  parameter  Pareto  distribution 
had  comparable  performance  to  the  WW  and  KK  distribution 
for  our  data  set. 
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